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induced  motion  noise  were  observed  in  the  spectra,  particularly  at  the  exposed  China  Lake 
sites.  Prominent  spikes  were  present  in  the  spectra  above  about  10  Hz;  except  for  several  of 
these  spikes  that  were  traced  to  the  tape  recorder  playback  electronics,  they  are  believed  to 
represent  genuine  man-made  ambient  signals.  Probable  sources  of  the  interference  include 
4-pole  synchronous  motor  emissions,  powerline  subharmonics,  and  electrified  railway  or 
telephone  ringing  signal  emissions.  In  California  the  coherence  between  the  nominally  parallel 
sensors  at  the  two  sites  was  relatively  high  (greater  than  0.75)  at  all  separations  and  showed 
definite  peaks  at  the  frequencies  of  the  Schumann  resonances.  The  coherence  was  highest  at 
frequencies  from  about  5  Hz  to  20  Hz,  usually  0.9  or  more.  No  structure  in  the  coherence 
attributable  to  the  effects  of  geological  inhomogeneities  was  seen  in  the  California  data.  The 
coherence  between  orthogonal  sensors  in  both  California  and  Norway  was  usually  0.25  or 
less  and  was  highly  varlablealn  Norway  coherence  at  9.3  km  approached  a  value  of  0.9  at 
the  highest  but  showed  stricture  that  probably  is  associated  with  the  conductivity  anomaly 
caused  by  a  nearby  fjord. 

An  algorithm  first  described  by  Fowler,  Smith,  a.id  Bostick  (FSB)  was  investigated 
analytically  for  its  ability  to  suppress  the  ambient  noise  by  using  frequency  domain  signal 
differencing.  The  FSB  processing  of  the  data  consistently  produced  a  difference  signal 
whose  noise  spectral  density  was  very  near  the  SQUID  magnetometer  intrinsic  noise  except 
for  a  noticeable  Increase  below  5  Hz.  Seismic-induced  motion  noise  in  the  SQUID  sensors 
probably  accounts  for  this  increase  in  the  processed  difference  below  5  Hz.  The  FSB 
processing  typically  produced  a  decrease  of  30  to  40  dB  or  more  compared  with  the  spectral 
amplitude  of  the  ambient  noise. 

The  results  reported  here  indicate  that  atmospheric  and  geomagnetic  noise  in  the  band 
from  0.3  to  40  Hz  can  be  reduced  essentially  to  the  instrument  noise  level  by  using  difference 
processing  with  two  magnetometers. 
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SPATIAL  COHERENCE  MEASUREMENTS  AND  EVALUATION  OF  A  NOISE 
REDUCTION  TECHNIQUE  FOR  AMBIENT  NOISE  FROM  0.3  TO  40  HZ 

SUMMARY 

Measurement*  of  the  ambient  magnetic  field  in  the  frequency  range  from  0.3  Hz  to 
40  Hz  have  been  made  at  a  mid-latitude  location  in  winter  (China  uake,  California)  and 
an  auroral  location  in  summer  (Tromsd,  Norway).  At  each  location  two  sensitive  triaxial 
SQUID  magnetometers  were  installed  at  sites  whose  separation  varied  from  0,5  km  to  15  km 
(in  California)  and  from  0.5  km  to  9.3  km  (in  Norway).  Approximately  two  weekB  of  nearly 
continuous  data  were  recorded  at  each  location.  The  data  analysis  in  this  report  is  concerned 
primarily  with  frequency  domain  investigations:  spectral  features  in  the  recorded  magne¬ 
tometer  signals  are  identified,  the  frequency-dependent  spatial  coherence  is  calculated  for 
various  combinations  of  sensor  axes  and  its  features  analyzed,  and  a  frequency  domain 
signal  differencing  technique  is  assessed  for  its  ability  to  suppress  the  ambient  noise.  The 
major  effort  was  expended  in  analyzing  and  assessing  the  noise  suppression  technique.  Our 
specific  comments  and  conclusions  regarding  each  area  of  the  data  analysis  are  as  follows: 

•  Ambient  noise  spectra.  Schumann  resonances  with  a  particularly  high  ratio  of  signal- 
to-instrument  noise  were  observed.  Six  orders  of  the  resonances  were  typically  observed  in 
the  Bpectra.  Significant  peaks  attributed  to  wind-induced  motion  noise  were  observed  in  the 
spectra  particularly  at  the  exposed  China  Lake  sites.  Prominent  spikes  were  present  in  the 
spectra  above  about  10  Hz;  except  for  several  of  these  spikes  that  were  traced  to  the  tape 
recorder  playback  electronics,  they  are  believed  to  represent  genuine  man-made  ambient 
signals.  Sources  of  the  interference  include  4-pole  synchronous  motor  emissions,  powerline 
subharmonics,  and  electrified  railway  or  telephone  ringing  signal  emissions. 

•  Spatial  Coherence.  In  California  the  coherence  between  the  nominally  parallel 
sensors  at  the  two  sites  was  relatively  high  (greater  than  0.75)  at  all  separations  and  showed 
definite  peaks  at  the  frequencies  of  the  Schumann  resonances.  The  coherence  was  highest  at 
frequencies  from  about  6  Hz  to  20  Hz,  usually  0.9  or  more.  No  structure  in  the  coherence 
attributable  to  the  effects  of  geological  inhomogeneities  was  seen  in  the  California  data. 
Larger  separation  distances  than  15  km  would  be  necessary  to  observe  geological  effects. 

A  large  decrease  in  the  coherence  at  certain  frequencies  was  noticeable  for  data  contami¬ 
nated  by  wind-induced  motion  noise,  The  coherence  between  orthogonal  sensors  in  both 
California  and  Norway  was  usually  0,25  or  less  and  was  highly  variable.  The  Norway 
coherence  results  revealed  a  surprisingly  low  coherence  of  about  0.5  between  sensors 
assumed  to  be  parallel  at  the  0,5  km  separation.  This  low  coherence,  taken  together  with  the 
difference  processing  results,  can  be  explained  in  a  consistent  manner  by  postulating  that 
the  remote  site  magnetometer  was  inadvertently  installed  with  the  horizontal  axes  rotated 
by  30°.  The  coherence  at  9.3  km  approached  a  value  of  0,9  at  the  highest  but  showed 
structure  that  we  believe  is  associated  with  the  conductivity  anomaly  caused  by  a  nearby 
fjord. 
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•  Fowler-Smith-Bostick  (FSB)  Processing.  The  FSB  algorithm  was  first  investigated 
analytically  and  an  expression  derived  from  the  processed  difference  signal  that  aided  in  the 
interpretation  of  the  processing.  The  FSB  processing  of  the  data  consistently  produced  a 
difference  signal  whose  noise  spectral  density  was  very  near  the  SQUID  magnetometer 
intrinsic  noise  except  for  a  noticeable  increase  below  5  Hz.  Seismic-induced  motion  noise  in 
the  SQUID  sensors  probably  accounts  for  this  increase  in  the  processed  difference  below 
5  Hz.  The  FSB  processing  typically  produced  a  decrease  of  10  to  30  dB  in  the  difference 
spectral  amplitude,  compared  with  the  simple  first  difference,  and  produced  a  decrease  of 
40  dB  or  more  compared  with  the  spectral  amplitude  of  the  ambient  noise. 

The  results  reported  here  indicate  that  atmospheric  and  geomagnet  c  noise  in  the  band 
from  0.3  to  40  Hz  can  be  reduced  essentially  to  the  instrument  noise  level  by  UBing 
difference  processing  with  two  magnetometers. 


INTRODUCTION 

This  report  is  concerned  with  the  measurement  and  characterization  of  naturally 
occurring  ambient  electromagnetic  noise  in  the  frequency  range  from  0.3  Hz  to  40  Hz  and 
with  the  evaluation  of  a  technique  to  reduce  the  noise  spectral  density.  The  lower  decade  of 
this  frequency  range  lies  in  the  ultra  low  frequency  (ULF,  0  to  3  Hz)  band;  noise  in  this 
band  is  generated  by  processes  that  occur  in  the  magnetosphere  and  upper  ionosphere.  At 
frequencies  above  3  Hz  in  the  extremely  low  frequency  (ELF,  3  to  3000  Hz)  band,  the 
noise  is  dominated  by  atmospheric  impulses  generated  by  near  and  distant  lightning  strokes. 
The  study  of  noise  in  this  frequency  band  and  techniques  for  its  mitigation  is  important  for 
submarine  magnetic  anomaly  detection  (MAD)  applications  and  for  proposed  ULF/ELF  low 
data  rate  communication  systems.  Noise  measurements  are  also  of  fundamental  interest 
as  a  ground-based  ionospheric/magnetospheric  diagnostic  technique  an  1  as  a  technique  for 
geophysical  exploration. 

Approximately  four  weeks  of  nearly  continuous  recordings  were  taken  at  sites  in 
California  and  Norway.  At  each  site  two  triaxial  superconducting  quantum  interference 
device  (SQUID)  magnetometers  were  used  that  were  separated  by  distances  varying  between 
0.B  km  and  IB  km.  Because  of  the  very  large  quantity  of  data  collected,  the  emphasis  in  this 
report  is  on  the  analysis  and  interpretation  of  selected  intervals  of  these  recordings  to 
determine  frequency  spectra  and  the  frequency  dependent  Bpatial  coherence  function,  and 
to  evaluate  a  frequency  domain  noise  cancellation  technique;  in  short,  frequency  domain 
analysis  is  favored.  A  future  report  Is  planned  that  will  examine  the  temporal  characteristics 
of  the  recorded  noise. 

Although  much  observational  data  of  atmospheric  and  ionospheric/magnetospheric 
noise  exist  in  the  litorature  (see  Ref.  1  for  a  list  of  recent  references),  most  of  these  data  were 
taken  with  conventional  loop-type  antennas  whose  self-noise  level  has  been  demonstrated 
[2]  to  exceed  often  the  average  ambient  noise  level  in  the  0,3  to  40  Hz  band,  particularly  in 
the  2  to  8  Hz  portion  of  this  band.  Measurements  with  sensitive  SQUID  magnetometers  have 
been  reported  by  Fraser-Smith  and  Buxton  [1,2J  and  Zimmerman  and  Campbel’  [3] ,  but 
these  measurements  in  each  case  were  for  one  fixed  site  in  one  geographic  location.  The 
measurements  reported  here  are  the  first  of  ambient  noise  and  its  spatial  coherence  taken 
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with  SQUID  magnetometers  at  dual  sites,  The  coherence  enhancement  technique  that  we 
have  chosen  to  use  was  first  described  by  Fowler,  Smith,  and  Bostick  [4] ,  who  applied 
the  method  to  coil  magnetometer  data  gathered  by  them  in  a  frequency  range  from  10“  4 
to  0.1  Hz.  However,  the  remote  reference  magnetotelluric  technique  of  GainM®,  et  nl 
[5,6]  has  some  resemblance  to  the  processing  techniques  described  in  •  his  report. 

In  the  next  section  the  various  factors  that  affect  the  spatial  cohere  ice  of  0.3  to  10  Hz 
noise  are  discussed,  and  an  algorithm  by  which  the  apparent  coherence  can  be  enhanced  is 
derived.  The  relationship  of  the  algorithm  to  linear  systems  theory  is  demonstrated  in  order 
to  assist  in  the  interpretation  of  the  processing  results.  Then  details  of  the  data  collection 
effort  are  given.  The  next  section  discusses  the  various  data  analysis  techniques  for  deter¬ 
mining  the  spectra  and  the  spatial  coherences,  and  for  investigating  the  coherer.ce  enhance¬ 
ment  algorithm.  An  expression  for  the  noise  spectral  density  after  the  coherence  enhance¬ 
ment  technique  is  applied  is  derived.  Next  the  results  of  the  data  analysis  are  given,  and 
finally  the  conclusions  of  the  study  and  recommendations  for  further  work  are  reported. 


SPATIAL  CORRELATION  AND  NOISE  REDUCTION 
Introduction 

Since  the  free  space  wavelengths  in  the  0.3  to  40  Hz  range  vary  between  10:1  Mm  and 
7.5  Mm,  one  might  expect  the  noise  fields  to  be  highly  correlated  over  very  large  distances. 
However,  on  the  earth's  surface  high  correlation  exists  over  distances  of  only  tenB  of  kilo¬ 
meters,  for  the  following  reasons: 

•  The  measurement  sites,  particularly  for  geomagnetic  noise  below  several  hertz,  are 
necessarily  in  the  complex  near  field  of  the  sources; 

•  The  sources  can  be  dispersed,  so  that  the  wavefronts  are  not  perfectly  coherent; 

•  Inhomogeneous  geology  in  the  vicinity  of  the  measurement  sites  can  distort  (or, 
from  a  different  point  of  view,  scatter)  the  incident  noise  fields,  thereby  lowering  the 
correlation. 

These  decorrelating  factors  are  all  intrinsic  properties  of  the  noi;  «  fields  themselves;  of  the 
three  the  last  factor  probably  contributes  the  most  to  the  lowi  mg  of  the  coherence.  Two 
additional  factors  associated  with  the  detection  of  the  fields  also  cause  an  apparent  lowering 
of  the  correlation: 

•  Misalignment  or  uncompensated  motion  of  the  presumed  parallel  vector  sensors  will 
lower  the  correlation; 

•  Intrinsic  sensor  noise  can  appreciably  degrade  the  measured  correlation  if  the  signal 
to  noise  ratio  is  low. 

The  goal  of  a  noise  reduction  technique  is  to  enhance  the  apparent  site-to-site  spatial 
correlation  by  compensating  in  some  manner  for  the  decorrelating  factors  listed  above,  Of 
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more  fundamental  interest  is  the  determination  of  the  degree  to  which  each  of  the  listed 
factors  contributes  to  the  measured  coherence.  For  example,  separation  of  source  and 
instrumental  effects  from  geology  effects  would  be  important  if  spatial  coherence  surveys 
were  to  be  applied  to  geophysical  exploration.  We  have  found,  however,  that  separation 
of  the  factors  is  difficult  to  accomplish  in  measured  data,  a  point  to  be  discussed  in  detail 
later. 

We  now  describe  a  technique  first  derived  by  Fowler,  Smith,  and  Bostick  [4]  for 
enhancing  the  correlation,  This  technique  was  applied  to  all  of  our  data  and  was  extensively 
evaluated. 


Fowler,  Smith,  and  Bostick  (FSB)  Technique 
Background 

The  FSB  algorithm  begins  by  forming  a  linear  relationship  (transfer  function) 
between  the  three  orthogonal  magnetic  field  vector  components  measured  at  one  receiver 
site  (denoted  the  base  site)  and  a  single  vector  component  at  the  other  site  (denoted  the 
remote  site)  that  is  given  by 


A (f)  -  R(f)  -  A(f)X(f)  -  B(f)Y(f)  -  C(f)Z{f ),  (1) 


where  R  is  the  Fourier  transform  (FT)  of  the  sensor  output  of  the  remote  site;  X  is  the 
FT  of  one  horizontal  component  of  the  base  site;  Y  is  the  FT  of  the  other  orthogonal  hori¬ 
zontal  component  of  the  base  site;  Z  is  the  FT  of  the  vertical  component  of  the  oase  site; 

A  is  the  difference  signal;  A ,  B,  and  C  are  coefficients  to  be  determined;  and  f  is  the  fre¬ 
quency.  The  vector  component  R  in  general  can  be  in  any  fixed  arbitrary  direction  at  the 
remote  site,  i.e.,  it  is  not  assumed  that  the  direction  of  R  is  in  the  same  compass  direction 
as  X,  Y,  or  Z  at  the  base  site.  However,  during  the  experimental  measurements  we  usually 
attempted  to  align  the  three  orthogonal  sensors  of  the  remote  site  SQUID  magnetometer  so 
that  they  were  approximately  parallel  with  the  base  site  sensors.  When  we  want  to  refer  to 
one  of  these  remote  site  sensors,  we  will  use  the  notation  Rx ,  Ry ,  and  Rt ,  where  the  sub¬ 
script  denotes  the  axis  at  the  base  site  that  the  remote  site  is  nominally  parallel  to.  See 
Fig.  1. 

The  coefficients  A ,  B,  and  C  express  a  linear  frequency  dependent  relationship  between 
the  base  and  remote  sites;  they  form  what  FSB  call  a  "geology  filter"  in  reference  to  the 
assumption  that  local  inhomogeneities  in  the  subsurface  geology  are  the  dominant  contri¬ 
bution  to  the  particular  values  of  the  coefficients.*  In  a  MAD  application  the  coefficients 
would  be  evaluated  in  the  known  absence  of  a  perturbation  signal  at  the  remote  site  (the 
remote  site  is  assumed  to  be  the  sensor  attempting  to  detect  a  target).  The  coefficients  thus 
determined  would  then  be  used  to  form  the  difference  according  to  Eq.  (1)  at  a  later  time; 


*Tho  use  of  three  ieparete  letteri  for  these  coefficients  Instead  of  a  single  subscripted  variable  leads  to  some 
awkwardness  In  the  resulting  equations.  However,  we  choose  to  adopt  the  notation  of  FSB. 
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Fig,  1  —  Geometry  Mid  notation*  for  bu«  and 
remota  *it««,  Th»  *«nior  Rt,  although  recorded 
part  of  the  time,  wa*  not  uaed  In  any  of  the 
data  proeexiing  reported  here. 


because  the  value  of  A(f)  is  now,  we  expect,  much  less  than  R(f),  the  anomaly  signal  would 
be  enhanced  (the  implicit  assumption  is  made  that  the  base  site  is  sufficiently  far  removed 
so  that  it  does  not  detect  the  anomaly).  In  a  geophysical  exploration  technique  the  specific 
structure  otA(f),  *(/),  and  C(f)  would  be  used  to  gain  knowledge  of  the  subsurface  geology. 
The  effectiveness  of  the  FSB  approach  depends  on  the  robustness  of  Eq.  (1),  i-e.,  the  degree 
to  which  it  actually  expresses  the  true  relationship  between  the  measured  fields  at  the  base 
and  reference  sites.  Also,  to  be  effective  as  a  MAD  technique,  the  values  of  A,  B,  and  C  must 
be  time-invariant;  by  scanning  the  list  of  decorrelating  factors  previously  given,  we  see  that 
time  invariance  requires  that  the  source  near-field  structure  must  have  a  small  effect  on  the 
coefficients  (otherwise,  the  coefficients  determined  for  one  source  configuration  would 
possibly  not  work  for  a  later,  different  configuration),  the  geological  structure  must  be 
essentially  invariant  (probably  an  excellent  assuinption),  and  the  BenBors  must  be  stationary. 

Coefficient  Evaluation 

To  estimate  the  coefficients  A(f),  B(f),  and  C(f)  in  Eq.  (1),  we  assume  that  the  input 
time  series  B(f),  X(t),  Y(t).  and  Z(t)  are  digitized,  and  a  total  of  N  sequential  discrete 
Fourier  transforms  are  computed.  These  N  FTs  may  be  computed  with  some  overlap  in  the 
input  time  series.  At  each  discrete  frequency  fa  within  the  FT  domain,  we  examine  '//,  the 
squared  difference: 


I  n  l 

-  L  R((f)~ Mnxtf)- mnYfif)- cv)z.(f)\ 

i  - 1 L  -* 

X  R*(f)  -  A*(f)X*(f) - B\f)Y*(f ) -cV)Z,*(flj  . 
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The  squared  difference  is  clearly  the  logical  error  quantity  for  the  determination  of  the 
coefficients,  In  the  usual  manner  then  we  set 


d\p 

d\p 

9R  *{A(f)) 

'  aim (A(f)} 

b\j> 

bit 

blm{B(f)) 

d\p 

3Rfl{C(f)} 

aim  {C(f)} 

where  Re  {•}  and  Im  {•}  denote  the  real  and  imaginary  parts,  respectively.  Applying  the 
operations  of  Gq.  (3)  on  the  expression  for  A  given  in  Gq.  (1)  yields  (deleting  the 
frequency  dependence) 


IM(X*  -  AZXtX*  +  BSYtX*  +  CXZ,X *  ( 

2/?, Y*  “AZXl Y*  +  Bl Y. Y*  +  CX,Zi Y*  ,  (4 ) 

2 RtZ*  •  A2X.Z*  +  BlYjZ*+  Clztf  . 


The  summations  involve  auto-  and  cross-power  density  spectra  between  all  combinations 
of  the  four  input  time  series.  The  solution  to  Eq.  (4)  can  be  written  in  matrix  form  (a  bar 
ib  used  to  denote  the  summation)  as 


at  each  discrete  frequency  in  the  FT  domain.  Equation  (5)  is  the  solution  we  seek  for  the 
coefficients  A,  B,  and  C  in  terms  of  quantities  that  can  be  evaluated  from  the  four  input 
time  series.  By  defining 
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r  “  m  b  cf , 


p 


(6) 


where  the  T  indicates  the  transpose.  Equation  (6)  can  be  written  compactly  as 


T  “  M“ 1  /?  . 


(7) 


The  explicit  equations  for  the  coefficients  A ,  B,  and  C  can  be  written  out  from  Eq.  (5) 
using  Cramer’s  rule  [7] .  The  resulting  expressions  are  very  involved  in  the  general  three- 
dimensional  case  ( X ,  y,  Z  0).  However,  the  essential  physics  of  the  FSB  noise  reduction 
technique  can  bo  examined  by  considering  only  the  two  horizontal  axes  at  the  base  site. 
Setting  Z  »  0  in  Eq.  (6)  and  solving  for  A  and  fl  gives  the  following  equations: 


aR  [coh(f?X)  -  coh (R  y)coh(  yxj] 
ax  |l  -  |coh(*y)|2] 

and 

[coh(fly)-  coh(i?^)coh(yX)J 

B  Oy  [l  -  |coh(Xr)|2] 

where  we  have  introduced  the  coherence  function  defined  by 


coh(cb) 


(<iP5F*)1/2 


(8) 


(9) 


(10) 


Hr-'- 
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The  identifications 


for  the  one-sided  average  spectral  noise  densities  have  also  been  made  in  Eqs.  (8)  and  (9). 
The  coherence  function,  as  we  have  defined  it  here,  is  a  complex  variable  whose  magnitude 
lies  between  0  and  1.0  and  is  an  indication  of  the  phase  correlation  between  the  two  signals. 
A  value  of  1.0  for  |coh(ab)|  indicates  perfect  correlation,  while  a  value  of  0.0  represents  no 
correlation. 

Equations  (8)  and  (9)  reveal  that  problems  will  occur  when  the  X  and  V  base  site 
signals  are  highly  correlated  with  each  other,  since  A  and  B  become  infinite  when 
coh(xy)  ■  1.  In  the  data  the  base  site  cross  coherences  were  generally  found  to  be  low, 
except  for  Isolated  tonal  interference  associated  with  man-made  signals  (such  as  powerline 
frequencies), 

Alternative  Derivation  of  the  Coefficient t 

The  above  derivation  of  Eqs.  (5),  (6),  and  (7)  for  evaluating  A ,  Bt  and  C  is  exactly  as 
given  by  FSB  in  [4] .  However,  by  viewing  the  determination  of  the  coefficients  in  a  some¬ 
what  different  manner,  the  equations  are  seen  to  be  identical  to  the  equations  commonly 
UBed  in  evaluating  linear  system  transfer  functions  (as  given,  for  example,  in  Bendat  and 
Piersol  [8] ).  The  demonstration  of  this  equivalence  is  important,  because  equations  for  the 
statistical  error  for  the  transfer  function  are  available  in  [8]  and  can  be  used  directly  to 
estimate  errors  in  the  determination  of  A ,  J3,  and  C.  These  statistical  error  equations  assist 
in  the  interpretation  of  the  spectral  features  of  the  coefficients  by  allowing  statistically 
significant  variations  to  be  identified. 

Equation  (1)  can  be  written  as 


R  -  A  »  AX  +ST  +  CZ, 


(ID 


which,  when  A  is  forced  to  zero,  is 
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R-AX  +  BY  +  CZ.  (12) 


Hence,  R  appears  as  the  output  of  a  linear  system  whose  inputs  are  X,  Y,  and  Z,  and  whose 
transfer  functions  are  A,  B,  and  C  respectively,  Reference  8  discusses  the  manner  in  which 
the  transfer  functions  of  a  multiple-input  linear  system  are  determined  and,  i’  irticular, 

emphasizes  the  precautions  that  must  be  observed  when  coherent  inputs  are  ,  wit.  To  ) 

demonstrate  the  connection  between  the  FSB  algorithm  and  the  equations  of  18] ,  we  will  \ 

for  simplicity  set  2  “  0  so  that  we  consider  only  a  two  "input”  system  ( X  and  Y).  From  [8]  J 

the  equations  for  A(f)  and  B(f)  can  be  written  js 


f 

l 


A(f) * 


SRx.y(f) 


sx x.yif)  ' 


and 


B(f)~ 


3f?  y.x 


(f) 


SYY.xW 


in  which  tiie  Sjj.k  are  termed  residual  spectra  and  are  defined  by 


SRl-J=‘SRJ 


and 


W-X’Y), 


(13) 


(14) 


(15) 


(16) 


where 

7%y  -  lcoh(XF)|2  . 

The  notation  has  been  taken  over  verbatim  from  (8] ;  the  Sj ,  are  the  single  sided  auto- 
and  cross-power  spectral  densities.  However,  note  that  in  terms  of  the  notation  we  have  used 
in  Eq.  (5), 

j 

Sjj  —  OjOj  coh(/«/)  {I,J-X,Yjl),  (17) 

'i 
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and 


Sjj  =  a\  (I  =  X,Y) .  (18) 

When  Eqs.  (17)  and  (18)  are  substituted  into  Eqs.  (13)  and  (14)  and  rearranged,  Eqs.  (8) 
and  (9)  (the  FSB  equations)  result.  Hence,  the  equivalence  between  the  FSB  method  and 
the  evaluation  of  linear  system  transfer  functions  has  been  demonstrated.  This  equivalence 
will  be  exploited  below. 

Expected  Noiie  Level  of  the  Difference  Signal  A 

The  expected  (theoretical)  noise  level  of  the  processed  difference  signal  A  1b  of  interest 
both  for  demonstration  of  the  noise  reduction  capabilities  of  the  FSB  algorithm  and  for 
comparison  with  the  experimental  results.  In  Appendix  A  an  expression  is  derived  for  the 
noise  power  in  the  processed  difference  signal,  denoted  by  o\ ,  for  the  two-dimensional  case. 
ThiB  equation  is  given  by  Eq.  ( A-4),  which  we  repeat  here: 

„  „  A  "  '•'RX  "  lRY~  Ixy  “  2Re]T/}x7/eyT'xyl 

-  -  .  (19) 

'EY 

where 


yjj  -  [coh(W)] 2  (I m  R,X\J  m  X,Y) . 

Equation  (19)  indicates  the  power  in  the  processed  difference  signal  is  equal  to  the  power 
in  the  reference  site  sensor  signal  modified  by  the  factor  in  brackets.  Three  special  cases  can 
be  examined  that  show  the  features  of  Eq,  (19). 

1.  yxy  1.0.  In  this  case  o\  grows  without  limit,  again  emphasizing  the  instability  in 
the  processing  when  highly  coherent  signals  are  present  in  both  base  site  axes.  More 
typically  7xy  0. 

2.  7 rx  —  I.Try  —  0, 7 xy  “  0.  This  case  corresponds  to  a  remote  sensor  that  is 
parallel  to  the  base  site  X-axis  in  the  presence  of  small  decorrelating  effects.  Then  Eq.  (19) 
becomes 

°A  -***)•  <2°) 

Thus,  the  higher  is  the  R  and  X  coherence,  the  lower  is  the  noise  in  the  difference  signal. 
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3.  yn x  “  Jr  y  “  ® .  7x  Y  m  **>y  value.  Thu.  case  corresponds  to  an  attempt  to  use 
completely  uncorrelated  signals  at  the  base  and  reference  sites  to  reduce  the  power  in  A.  As 
expected,  in  this  case,  i.e.  no  noise  reduction  occurs. 

Appendix  A  carr.es  Eq.  (19)  further  in  order  to  obtain  an  expression  for  in  terms 
of  both  the  SQUID  sf  nsor  noise  level  (which  one  intuitively  expects  to  be  the  ultimate 
limit  for  o^)  and  contributions  to  the  noise  that  might  not  be  cancelled  by  the  construction 
of  the  algorithm  (arising,  for  example,  from  nonlinear  effects  and  possible  uncorrelated 
noise  sources  other  then  the  intrinsic  sensor  noise).  The  simple  resulting  equation  is  given  by 

ffA  = 2on  + caR  •  (21) 


where  is  the  SQUID  sensor  noise  level  <  c  is  a  measure  of  the  fractional  deviation  of 
the  system  transfer  function  from  a  linear  response.  The  quantity  e  is  defined  with  greater 
care  in  Appendix  A.  Equation  (21)  indicates  that  the  power  in  the  difference  signal  is  within 
a  factor  of  two  of  the  SQUID  sensor  noise  level  (for  e  ■  0). 


DATA  COLLECTION 

The  data  discussed  in  this  report  were  collected  in  two  widely  separated  geographic 
areas.  Data  sample  periods  of  approximately  two  weeks  duration  were  obtained  during 
January  1978  at  the  Naval  Weapons  Center  at  China  Lake,  California  (Lat.  36°N,  Long. 
118°W)  and  during  June  197-8  at  Troms0,  Norway  (Lat.  68°N,  Long.  19°E).  The  China 
Lake  location  (Fig.  2)  consists  of  flat  and  sandy  desert  terrain  while  the  Norway  location 
(Fig.  3)  is  a  mountainous  fjord  region  within  the  auroral  zone.  The  base  and  remote  site 
separation  distances  used  at  each  area  were: 

•  California— 0.5,  4,8  and  14.6  km 

•  Norway— 0.5,  3.5  and  9.3  km 

The  receiving  systems  employed  at  each  site  were  nearly  identical.  The  magnetometers 
used  superconducting  quantum  interference  devices  (SQUIDs)  to  achieve  a  low  intrinsic 
noise  level  and  wide  dynamic  range.  The  SQUIDs  were  cooled  by  liquid  helium.  Table  1  lists 
the  specifications  of  the  four  systems  used  in  these  measurements  and  indicates  the  site  at 
which  each  was  used.  The  principles  of  operation  of  a  SQUID  will  not  be  discussed  in  this 
report,  but  the  reader  who  is  interested  in  more  details  is  directed  to  [9] . 

The  triaxial  SQUIDs  were  sensitive  to  three  orthogonal  directions— two  horizontal  axes 
and  a  vertical  axis.  The  SQUIDs  were  installed  with  the  two  horizontal  directions  along 
mugnetic  north-south  and  east-west  axes.  A  compass  and  spirit  level  were  used  to  perform 
the  azimuthal  alignment  and  establish  the  vertical  axis.  The  SQUID  magnetometers  were 
buried  in  the  soil  as  deeply  as  possible  to  anchor  them  firmly  and  to  reduce  their  wind  cross 
section.  A  low  wooden  shelter  was  used  to  shield  the  magnetometers  from  wind  and  the 
weather.  Figure  4  is  a. photograph  of  the  installation  in  California. 
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Fig.  2  —  Location  of  recording  sites  at  China  Lake,  California 
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Table  1  —  SQUID  Magnetometer  Systems 


System 
(All  Dynamic 
Ranges^  130  dB) 

Site 

Approximate 
Intrinsic 
Noise,  nG/y/Hz 

Internal 

RFI 

Shielding'*' * 

Liquid 

helium 

fill 

Interval, 

days 

NRL-owned  ELF 
system  [10] 

Calif. 

base 

0.10 

Two  layers, 
f0  -  BO  Hz 

7 

SHE  Corp  leased 
system  I 

Calif. 

remote 

0.24 

Single  layer, 
fo  -  35  Hz 

7 

SHE  Corp  leased 
system  II 

Norway, 

base 

<0.4 

Single  layer, 
f0  “  37  Hz 

7 

System  borrowed 
from  NCSC* 

Norway, 

remote 

0.20 

_ 

None 

5 

*NCSC  »  Naval  Coastal  Systems  Center,  Panama  City,  Fla, 
**/D  ■  3dB  cornor  frequency  for  RFI  ahielda. 


Fig.  4  —  Typical  installation  of  SQUID  magnetom¬ 
eter.  A  wind  acrecn  Is  placed  over  the  top  or  the  mag¬ 
netometer  during  the  recording  periods, 
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In  both  Norway  and  California  it  was  our  intention  to  align  the  axes  at  the  separated 
sites  in  approximately  the  same  directions,  so  that  in  the  analysis  R  would  be  essentially 
parallel  to  X  or  Y  at  the  base  site  (depending  on  which  axis  at  the  remote  site  is  chosen).  It  is 
not  essential  for  the  coherence  enhancement  technique  that  R  be  parallel  to  one  of  the  base 
site  axes,  since  the  coefficients  in  Bq.  (5)  in  effect  perform  a  coordinate  rotation.  However, 
for  comparison  with  a  direct  difference  (denoted  the  first  difference)  and  for  direct  compu¬ 
tation  of  the  coherence  between  physically  parallel  axes  for  the  two  separated  sites,  we 
attempted  to  align  the  two  sets  of  sensors  as  parallel  as  possible.  (At  the  0.5  km  separation 
in  Norway,  however,  the  axes  were  inadvertently  misaligned  by  an  angle  of  about  30°.) 

The  signal  conditioning  and  data  recording  systems  are  shown  in  Fig.  5.  The  SQUID 
field  measurement  noise  level  was  preserved,  and  the  signal  to  noise  ratio  (SNR)  established, 
by  the  low  noise  Princeton  Applied  Research  Model  113  preamplifier;  this  SNR  was  not 
degraded  by  the  remainder  of  the  signal  path.  Following  the  preamplifier  were  power 
frequency  notch  filters  that  used  modules  manufactured  by  Burr  Brown,  Model  ATF76 
"cries  circuits,  with  Q  “  10.  The  configuration  consisted  of  three  cascaded  filters,  with  one 
at  the  primary  power  frequency  (50  or  60  Hz)  and  the  other  two  staggered  by  1  Hz  above 
and  below  the  primary  power  frequency,  In  addition,  a  single  filter  was  used  at  both  the 
first  and  second  frequency  harmonics,  The  amplified  and  filtered  signals  at  the  remote  site 


REMOTE  SITE 


Fig.  S  —  Signal  conditioning  and  recording  electronic* 
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were  input  directly  to  an  analog  FM  UHF  telemetry  system  operating  at  a  carrier  frequency 
of  229.9  MHz.  The  telemetered  remote  site  signals  and  the  base  site  signals  were  then  low 
pass  filtered  by  a  filter  with  a  two  pole  Butterworth  response  with  a  40-Hz  comer 
frequency,  for  anti-aliasing  purposes,  before  being  recorded. 

Three  types  of  data  recording  were  used.  A  hard  copy  was  made  using  a  pressure  ink 
strip  chart  recorder,  and  magnetic  tape  recordings  were  produced  using  both  an  analog 
14-track  instrumentation  tape  recorder  and  a  12-bit  digital  tape  recorder. 


DATA  ANALYSIS  TECHNIQUES 

The  recorded  data  tapes  were  scanned,  and  intervals  that  appeared  to  be  free  of  any 
unusual  man-made  Interference  were  selected  for  analysis.  The  digital  data  tapes  recorded 
in  California  were  used  directly;  howevet,  since  no  digital  recorder  was  available  in  Norway, 
the  analog  tapes  recorded  there  were  digitized  in  the  laboratory  using  a  12  bit  analog-to- 
digitol  converter.  The  digitizing  rate  in  both  California  and  Norway  was  83.3  Hz  for  each 
magnetometer  axis,  corresponding  to  a  time  sample  interval  of  12  ms.  The  periods  selected 
for  analysis  are  listed  in  Table  2. 


Table  2  —  Data  Intervals  Analyzed  in  this  Report 


Location 

Separation  (km) 

Analysis 

Interval 

(Day/h/min) 

Analysis 
Results  in 
Figures 

Remarks 

China  Lake, 
California 

0.5 

12/21/40 

-12/21/50 

6,12,17,18 

4.8 

16/8/0 

-15/8/10 

7 

14.6 

18/23/20 

-18/23/30 

9,13,19 

Calm 

conditions 

20/13/40 

-20/13/50 

8,14,20 

High  winds 

Tromsp, 

Norway 

0.5 

173/12/20 

-173/12/30 

10,15,22 

9.3 

176/19/30 

-176/19/40 

11,16,  23 
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Spectral  Analysis 

We  used  a  standard  discrete  Fast  Fourier  transform  (FFT)  program  to  generate 
sequential  1024  real  point  spectra  for  each  of  the  magnetometer  time  series.  Each  spectrum 
required  a  total  of  (1024)  (.012)  *  12.29  s  of  data.  We  employed  a  Hanning  window  on  each 
12.29-s  time  slice,  which  we  chose  for  its  compromise  between  spectral  resolution  and 
leakage.  The  successive  time  slices  were  overlapped  by  50%,  so  that  half  of  the  data  in  each 
time  slice  was  “new"  and  half  was  “old".  This  60%  overlap  is  common  in  time  series  analysis 
and  has  been  shown  [11]  to  reduce  the  error  in  spectral  estimation.  The  spectra  were 
computed  over  a  bandwidth  of  0  to  40  Hz;  this  upper  frequency  is  just  below  the  Nyquist 
frequency  of  (83.33/2)  or  41.7  Hz. 

Although  the  frequency  response  of  the  SQUID  sensor  is  flat  to  at  least  1  kHz,  the 
eddy  current  electromagnetic  interference  (EMI)  shields  and  the  low  pass  filter  introduced 
an  amplitude  and  phase  distortion  that  had  to  be  compensated.  The  frequency  response  of 
each  SQUID  magnetometer  was  measured  using  an  external  drive  coil  to  generate  a  known 
magnetic  field  at  the  sensors.  This  measured  response  was  then  convolved  with  the  recorded 
data  in  order  to  compensate  for  the  shield  and  filter  high  frequency  roll-off. 

Confidence  limits  for  the  spectra  were  derived  in  the  following  manner.  The  effective 
spectral  bandwidth  Bt  is  given  by  Bt  -  (83.3  Hz/1024)(r/2)  -  0.128  Hz,  where  the  ir/2 
factor  is  necessary  to  obtain  the  effective  noise  bandwidth  [8] .  According  to  Nuttal  [11] , 
the  equivalent  degrees  of  freedom  Nt  for  60%  overlapped  spectra  for  a  total  record  length 
of  T  s  and  a  Hanning  window  is  given  by 


Nf  -2.71BaT-  3.36.  (22) 


Given  a  value  of  Nt ,  the  80%  confidence  interval  spread  r?  is  given  in  dB  by  [8] 


V 
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dB. 


(23) 


Combining  Equations  (22)  and  (23)  yields 
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or,  for  Be  -0.128  Hz, 
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V0.347T-  4.36 


(25) 


For  the  T  -  593  s  of  data  used  in  most  of  the  spectra  shown  below,  the  value  of  Nt  m  201, 
and  the  80%  confidence  spread  is  1.13  dB. 


Coherence  Analysis 

The  coherence  between  the  nominally  parallel  sixes  at  the  base  and  remote  sites  was 
calculated  from  the  equation 


H/)^(F|  _ 


(26) 


where  8(f)  and  R  (f)  are  the  FTs  of  the  nominally  parallel  a.,  ’s  at  the  base  and  remote  sites, 
respectively.  The  coherence  is  a  time  averaged  quantity  and  ;*  *  measure  of  the  time  stability 
of  the  relative  phase  between  S(f )  and  R(f),  The  cross-coherence  between  the  base  site  X 
and  Y  sensors  was  also  calculated  using  an  equation  similar  to  Eq.  (26).  The  average  for  the 
coherence  functions  presented  later  was  carried  out  over  593  s  of  data. 

Confidence  limits  for  the  coherence  can  be  derived  according  to  the  equations  given  in 
[8] .  For  an  estimated  value  of  coherence  7  2  and  Nt  equivalent  degrees  of  freedom,  the  true 
value  of  coherence  y2  falls  in  the  100(1  -  a)  percent  confidence  interval  defined  by 


tanh2  { - 


1 

Nt-  2 


(27) 


<  tanh 


where  Z0  is  the  50a  percentage  point  of  a  Qausslan  distribution,  and 


1  1+7 

{•  ■  -  In — ~r 
5  2  1-7 


(28) 
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The  same  equation  for  Ng  can  be  used  here  as  was  used  in  deriving  the  spectral  confidence 
limits,  since  the  coherence  plots  below  were  also  obtained  from  60%  overlapped  and 
windowed  spectra,  Table  3  lists  the  80%  confidence  limits  for  the  coherence  derived  for 
Be  *  0,128  Hz,  T  ■  593  s,  and  Ne  -  201.  The  confidence  interval  decreases  as  the  coherence 
increases. 


Table  3  —  Coherence  80%  Confidence  Limits 


Estimated 

0.2 

0.4 

0.6 

0.8 

0.9 

coherence  7  2 

Lower  limit 

0.13 

0.33 

0.64 

0.76 

0.88 

Upper  limit 

0.26 

0.46 

0.66 

0.83 

0.91 

FSB  Analysis 

The  FSB  algorithm  was  evaluated  using  the  60%  overlapped  spectra  as  the  input  data. 
Equation  (5)  was  implemented  on  the  Naval  Research  Laboratory  Advanced  Scientific 
Computer  (ASC)  in  a  straightforward  manner  to  evaluate  the  frequency  dependent  param¬ 
eters  A,B,  and  C  for  eaah  of  the  data  intervals  given  in  Table  2. 

Reference  8  discusses  both  the  bias  (systematic)  error  and  the  variability  (random) 
error  associated  with  measuring  system  transfer  functions.  Since  we  have  shown  that 
evaluation  of  the  coefficients  in  the  FSB  algorithm  is  equivalent  to  measuring  a  system 
transfer  function  (where  R  is  the  output  and  X ,  Y,  and  Z  are  the  Inputs),  we  can  use 
without  modification  the  equations  given  in  [8] .  The  error  equations  for  inputs  assumed  to 
be  uncorrelated  are  relatively  straightforward,  whereas  the  existence  of  correlation  between 
the  inputs  requires  complicated  expressions.  We  will  use  equations  that  assume  the  inputs 
are  uncorrelated,  although  it  will  be  seen  below  that  there  is  a  small  degree  of  correlation 
between  them.  The  correlation  between  inputs  is  sufficiently  low  that  the  uncorrelated- 
input  expressions  provide  a  reasonable  estimate. 

The  bias  error  for  A  is  given  by 


\Mf)\  -  \A(f)\ 


(29) 


where  A(f)  is  the  coefficient  estimate  and  °Xn  is  the  extraneous  noise  power  portion  of  the 
auto-power  spectrum  XX*.  Similar  expressions  hold  for  B(f)  and  C(f).  For  the  data  in  this 
report  the  extraneous  noise  power  results  from  the  intrinsic  sensor  noise  of  the  magnetom¬ 
eter;  as  will  be  shown  below,  the  ratio  o\jo\  is  typically  0.01,  implying  that  the  bias  error 
is  negligible, 
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The  variability  error  is  expressed  as  a  1  -  a  confidence  interval  for  the  true  coefficient 
amplitude  |A|  and  is  given  by  the  following  expression  (in  which  the  frequency  dependence 
is  suppressed): 


|A|-  n<  |A|<  \A\  ,  (30) 

where 


and  F[2,  N  -  2;  a]  **  100  a  percentage  point  of  an  F  distribution  with  2  and  Ne  -  2 
degrees  of  freedom, 

The  other  variables  in  these  equations  have  the  same  interpretation  as  given  earlier. 
Similar  equations  hold  for  the  B  and  C  coefficient  variability  errors. 

The  confidence  interval  is  indicated  on  each  plot  of  A,  B,  ur  C  in  the  next  section.  As 
an  example  of  a  confidence  interval  computation  consider  a  typical  case  for  X  nominally  par¬ 
allel  to  R:  7»  ~  0,9,  (o/{  /cjj )  -  1.0,  Ne  -  200,  and  ft  -  0.0B.  Since  F[2,198;0,06]  -  3.04, 
the  confidence  interval  is  0,065  for  |A|.  Far  this  same  example  typically  “  °-3  and 
oj  /<j$  ■  1.0;  the  confidence  interval  for  |B|  is  then  0.14.  For  ICI  typically  y%H  =*  0.3  and 
nj |  ■  (8)a  (see,  for  example,  Fig,  6  and  compare  Bx  and  Z).  The  confidence  interval 

in  this  case  for  ICI  is  0.50,  substantially  larger  than  for  A  and  fl.  Note  that  the  confidence 
interval  can  vary  with  frequency,  since  both  y2  and  the  ratio  of  the  spectral  densities 
depend  on  frequency. 


RESULTS 

Spectral  Analysis  Results 

The  amplitude  spectra  for  the  measured  data  are  shown  in  Figs.  8  to  11.  These  Bpectra 
are  typical  for  each  geographic  area  and  for  each  site  separation  distance.  The  spectra  with 
labels  X  and  Y  refer  to  the  two  horizontal  orthogonal  axes  of  the  trlaxial  NQUID  magnetom¬ 
eters  and  the  Z  label  refers  to  the  vertical  axis  (See  Fig.  1),  The  Base  label  corresponds  to 
the  magnetometer  site  for  which  all  three  Bensor  signals  were  Incorporated  into  the  geology 
filter,  i.e.,  X  Base  refers  to  the  term  X(f)  in  Eq.  (1),  etc,  The  Remote  label  corresponds  to 
the  magnetometer  site  for  which  a  single  sensor  signal  was  incorporated  into  the  geology  fil¬ 
ter,  i.e.,  either  X  Remote  or  Y  Remote  refers  to  the  term  R(f)  in  Eq.  (1).  The  Z  Remote 
spectra  are  not  displayed  because  the  vertical  axis  of  the  atmospheric  noise  usually  contained 
the  least  amount  of  information  (relative  to  the  horizontal  axes);  hence,  when  an  equipment 
outage  such  as  telemetry  channel  failure  or  tape  recorder  record  electronics  malfunction 
occurred,  Z  Remote  was  always  the  signal  sacrificed. 
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Each  of  the  plots  displays  spectra  that  have  been  averaged  for  593  s  (about  10  min) 
of  data  using  50%  overlap.  The  time  given  in  these  figures  is  Universal  (UT);  to  convert  to 
local  time  (LT)  for  the  California  data,  LT  -  UT  -  8  and  for  the  Norway  data,  LT  ■  UT  +  1. 

The  most  prominent  features  of  the  spectra  displayed  in  Figs.  6  to  11  are  the  Schumann 
resonances.  These  resonances  have  a  fundamental  frequency  in  the  vicinity  of  7.5  to  8  Hz  [  12  ] 
and  harmonics  at  frequencies  of  approximately  14,  20,  and  26  Hz.  Varying  degrees  of  spec¬ 
tral  features  are  visible  at  frequencies  below  the  Schumann  peaks.  Figures  8  and  9  show  sev¬ 
eral  well  defined  peaks  between  0.5  and  2.0  Hz  that  result  both  from  micropulsation  activity 
and  minromotion  of  the  antennas  in  the  earth’s  magnetic  field.  In  particular,  the  peaks  near 
2  Hz  in  Figs.  6  and  7  result  from  motion  noise  of  the  vector  sensors  in  the  earth’s  magnetic 
field  caused  by  wind.  Figure  9  is  a  second  set  of  spectra  for  the  14.6  km  separation  in 
California  taken  during  a  period  when  there  waB  essentially  no  wind.  The  peaks  near  2.0  Hz 
are  absent.  Although  the  SQUID  magnetometer  was  firmly  sunk  into  the  ground  and  pro¬ 
tected  by  a  windscreen,  some  motion  of  the  magnetometer  was  apparently  unavoidable  in  a 
strong  wind. 

A  number  of  narrowband  tones  are  evident  in  the  spectra  above  about  10  Hz.  The 
prominent  tones  at  10.2  Hz  and  (sometimes)  at  20.4  Hz  were  caused  by  the  tape  recorder 
electronics.  Another  tone  sometimes  visible  at  23.2  Hz  has  also  been  traced  to  the  tape 
recorder  electronics.  The  remaining  tones,  we  feel,  are  caused  by  various  types  of  genuine 
electromagnetic  interference  (genuine  in  the  sense  that  the  SQUID  sensor  is  responding  to 
incident  signals;  they  are  not  artifacts  due  to  subsequent  processing).  The  tone  at  about 
30.6  Hz,  particularly  evident  in  the  California  base  spectra,  is  Bimilar  to  tones  frequently 
observed  by  other  workers  that  were  attributed  by  Lokken  [13]  to  emissions  from  4-pole, 
60-Hz  synchronous  motors.  A  similar  tone  at  26.2  Hz  in  the  Norway  data  may  result  from 
60-Hz,  4-pole  synchronous  motors.  Other  subharmonics  or  near  subharmonics  of  the  local 
powerline  frequency  are  sometimes  visible  (for  example,  13.2  Hz  in  the  Norway  data)  that 
may  possibly  be  a  result  of  a  phenomenon  described  by  Madden  and  Thompson  [14] .  They 
observed  subharmonics  of  the  local  powerline  frequency,  which  they  attributed  to  nonlinear 
interactions  between  the  generators  on  a  power  grid, 

The  tone  at  16.7  Hz,  observable  in  most  of  the  Norway  data,  is  worthy  of  special 
mention.  This  spike  is  caused  either  by  the  telephone  ringing  Bignal  or  by  emission  from 
Norwegian  AC  electrified  railways.  Since  both  frequencies  are  nearly  the  same,  it  is  difficult 
to  separate  the  two  sources  without  further  measurements. 

Finally,  in  anticipation  of  the  correlation  anrlysis  in  the  next  section,  the  correspond¬ 
ing  base  and  remote  parallel  axes  are  observed  to  produce  spectra  with  very  similar  features; 
qualitatively  speaking,  they  “look  alike." 


Coherence  Results 

Coherence  plots  for  the  data  samples  are  presented  in  Figs.  12  through  16.  Each  of 
these  figures  shows  the  coherence  between  the  two  horizontal  nominally  parallel  sensors  at 
each  site  and  the  cross  coherence  between  the  two  horizontal  orthogonal  sensors  at  the  base 
site  and  at  the  remote  site.  The  figures  Bhow  the  magnitude  of  defined  in  Eq.  (26).  Con¬ 
fidence  intervals  calculated  from  Eq.  (27)  are  shown  for  a  coherence  value  indicated  by  the 
point  within  each  interval. 
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Fig.  13— Coherence  between  Indicated  teniori  for  date  recorded  et  14.6  km  separetlon  at  China  Lake 
California.  The  80%  confidence  Interval  le  ehown. 
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Fig.  14  —  Coherence  between  indicated  leniora  for  data  recorded  at  14,6  km  aeparatlon  at  China  Lake, 
California  during  period  of  high  wlnda.  The  80%  confidence  Interval  la  ahown. 
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California  Data 

For  a  separation  of  0.5  km  the  coherence  shown  in  Fig.  12  below  20  Hz  for  both  the  X 
and  Y  parallel  axes  is  relatively  high  (greater  than  0,75).  A  peaking  of  the  coherence  is 
noticeable  at  the  Schumann  resonance  frequencies,  and  deep  nulls  or  peaks  are  evident  at 
the  frequencies  of  the  tonal  interference  present  in  the  spectra  of  each  sensor,  depending  on 
the  relative  phase  stability  of  the  tones  at  each  sensor.  Above  20  Hz  the  coherence  of  the 
parallel  sensors  decreases,  presumably  reflecting  the  decrease  in  the  coherence  of  the  noise 
between  the  two  sites  at  the  shorter  wavelengths. 

At  the  increased  separation  of  14.6  km  the  overall  coherence  between  the  parallel 
sensors  decreases  noticeably,  but  not  drastically.  Dips  again  are  evident  at  some  of  the  fre¬ 
quencies  of  the  tonal  interference  observed  in  the  spectra,  and  the  coherence  decreases  with 
increasing  frequency  above  20  Hz,  In  the  coherence  plot  for  data  recorded  during  the  high 
winds  (Fig.  14),  definite  structure  is  visible  below  about  4  Hz.  Especially  noticeable  is  the 
dip  at  a  frequency  of  1.1  Hz  in  the  X  axis  coherence.  The  low  coherence  is  probably  a  result 
of  the  motion  of  the  antennas  in  the  earth's  magnetic  field  that  is  uncorrelated  at  the  two 
sites,  in  Fig,  13  the  coherence  between  X  base  and  X  remote  is  about  0.9  and  shows  rela¬ 
tively  little  structure;  however,  the  Y  base  -  Y  remote  coherence  is  substantially  lower.  We 
are  not  certain  why  the  east-west  direction  produces  a  lower  coherence. 

The  coherence  between  the  orthogonal  X  and  Y  base  Bensors  is  relatively  low  and  vari¬ 
able,  having  an  average  value  of  about  0.25.  Some  of  the  higher  variability  at  the  lower  values 
of  coherence  is  due  to  the  wider  confidence  intervals  at  low  coherence  (Table  3). 

Norway  Data 

Figure  15  displays  the  coherence  between  the  X  axes  at  the  0.5  km  separation  in 
Norway.  The  coherence  at  all  frequencies  is  substantially  less  than  the  coherence  at  0.5  km 
measured  in  California,  and  a  general  decrease  of  the  coherence  with  increasing  frequency 
can  be  seen;  a  peaking  of  the  coherence  at  the  Schumann  resonance  frequencies  is  also  evi¬ 
dent.  A  somewhat  lower  coherence  in  Norway  was  not  unexpected,  since  the  mountainous 
terrain  should  cause  more  distortion  of  the  incident  noise  fields.  The  magnitude  of  the 
decrease  indicated  by  comparing  Figs,  15  and  12,  however,  is  much  greater  than  we  expected, 
particularly  at  the  relatively  close  0.5  km  separation. 

Before  ascribing  the  decreased  coherence  in  Fig.  1 5  to  the  effects  of  the  mountainous 
terrain,  Fig.  16  should  be  examined.  This  figure  shows  the  coherence  between  the  X  and  Y 
sensors  for  the  9.3  km  separation  in  Norway,  and  for  the  X  axes  indicates  that  the  coherence 
is  higher  than  the  coherence  at  0.5  km.  We  believe  that  the  relatively  low  coherence  at  0.6  km 
in  Norway  (in  comparison  with  the  9.3  km  coherence)  indicates  an  inadvertent  misalignment 
of  about  30  degrees  between  the  two  magnetometers.  That  is,  the  coherence  at  O.b  km  would 
be  higher  if  the  sensor  axes  were  indeed  parallel.  ThiB  hypothesis  i«  supported  by  the  results 
of  the  FSB  processing  discussed  in  the  next  section. 

The  X  -  Rx  and  Y-  Ry  coherence  at  9,3  km  in  Fig.  16  show  two  distinct  variations. 
The  X  -  Rx  coherence  has  a  peak  near  1 0  Hz,  falling  below  this  peak  at  higher  and  lower 
frequencies  (similar  to  the  California  coherences),  whereas  the  Y  -  Ry  coherence  decreases 
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with  increasing  frequency.  A  possible  explanation  of  these  two  variations  of  the  coherences 
in  terms  of  the  local  geology  will  be  given  in  the  next  section  when  the  FSB  analysis  is  dis¬ 
cussed. 

As  with  the  California  data,  the  coherence  between  the  base  site  X  and  Y  axes  is  rela¬ 
tively  low  and  variable. 


FSB  Processing  Results 

Figures  17  to  22  inclusive  give  the  results  of  processing  the  recorded  data  according  to 
the  FSB  method  outlined  above.  Each  figure  contains  the  following  plots:  the  coefficients 
A,  B,  and  C;  the  simple  first  difference;  the  processed  difference;  and  the  improvement  in 
the  difference  due  to  the  FSB  processing.  The  first  difference  is  the  difference  between  the 
selected  remote  sits  sensor  (Rx  in  all  cases  except  Fig.  18)  and  the  sensor  at  the  base  site 
that  is  nominally  parallel  to  it  (X  in  all  cases  except  Fig.  18).  The  improvement  is  defined  by 


Improvement  *  20  logj  0 


_ First  difference 

Processed  difference 


(31) 


All  plots  have  a  logerithmic  frequency  scale,  and  the  complex  magnitude  of  the  indicated 
quantity  is  plotted.  All  quantities  were  generated  from  693  s  of  data,  corresponding  to  ap¬ 
proximately  200  degrees  of  freedom.  Confidence  intervals  calculated  from  Eq.  (30)  are 
shown  on  the  plots  for  4,  B,  and  C  in  each  case.  If  the  interval  haB  an  appreciable  variation 
with  frequency,  it  is  shown  at  spot  frequencies  either  above  or  below  the  curve;  otherwise, 
one  interval  for  the  entire  plot  is  shown.  The  confidence  interval  assists  in  distinguishing 
between  significant  features  and  statistical  fluctuations  in  the  coefficient  spectral  variation. 

California  Data 

Figures  17  and  18  show  the  results  for  the  0.6  km  separation  in  California.  In  Fig.  17 
Rx  is  used  for  the  remote  site  axis,  and  in  Fig  18  Ry  is  used,  Note  that  in  these  plots,  and 
in  other  figures  below,  the  vertical  scale  factor  varies,  depending  on  the  maximum  value 
plotted,  so  that  casual  visual  comparison  of  the  coefficients  may  be  misleading. 

As  anticipated  for  this  relatively  short  separation,  in  Fig.  17  the  coefficient  A  (in  this 
case  the  transfer  function  between  the  parallel  sensors  at  the  two  sites)  is  very  close  to  unity 
and  displays  essentially  no  structure  or  variation  with  frequency.  The  average  value  of  A, 
shown  by  the  dotted  line,  is  0,87.  The  fact  that  the  variation  of  A  with  frequency  is  flat  im¬ 
plies  that  the  non-unity  value  of  A  is  probably  caused  by  a  uniform  gain  difference  between 
the  X  and  Rx  sensors,  rather  than  by  geological  homogeneities  (differences  related  to  geo¬ 
logical  inhomogeneities  would  probably  not  produce  a  uniform  effect  across  the  frequency 
band).  A  misalignment  error  can  be  ruled  out  by  the  high  coherence  in  Fig.  12,  since  an 
alignment  error  (unlike  a  gain  difference)  lowers  the  coherence.  Coefficient  B  is  low  in  am¬ 
plitude,  indicating  that  relatively  little  transverse  component  of  the  magnetic  field  at  the 
base  site  has  to  be  admixed  to  the  difference  between  the  nominally  parallel  magnetic  north- 
oriented  sensors.  The  spikes  in  both  coefficient  B  and  C  are  due  to  the  various  interference 
sources  identified  in  the  spectra  of  each  input. 
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The  processed  difference  in  Fig.  17  for  the  0.5  km  separation  is  very  near  the  average 
sensor  noise  level  of  2  X  10" 1 0  G/^/Hz  across  the  entire  frequency  band,  although  there  is 
a  perceptible  increase  in  the  processed  difference  signal  at  frequencies  below  3  Hz.  This 
increase,  as  will  be  seen  in  other  cases,  is  present  to  various  degrees  in  all  of  the  processed 
differences  and  represents  the  effect  of  a  site-to-site  decorrelation  factor  that  cannot  be  cor¬ 
rected  by  the  FSB  processing.  In  other  words,  it  arises  from  a  non  zero  value  of  e  in 
Eq,  (21).  Later,  we  shall  adduce  that  this  increase  is  possibly  the  result  of  uncorrelated 
motion-induced  signals  caused  by  microseismic  noise  at  the  two  sites.  The  improvement  in 
the  difference  signal  due  to  the  processing  ;*  plotted  in  Fig.  17  and  has  an  average  value 
across  the  band  of  about  9  dB.  A  comparison  with  the  raw  spectra  of  this  data  (Fig.  6) 
shows  that  the  processing  suppresses  ambient  noise  in  this  frequency  band  by  up  to  40  dB. 

Figure  18  shows  the  results  for  the  0.5  km  separation  in  California  when  Ry  is  taken  as 
the  single  sensor  at  the  remote  site.  As  expected,  the  coefficient  B  is  now  of  the  order  of 
unity  and  A  has  a  low  value.  The  average  value  of  B  shown  by  the  dotted  line  is  0.98;  this 
value  implies  that  any  gain  errors  in  the  time  series  data  are  very  small  for  the  A., and  Y  axes. 
As  in  the  case  of  the  processing  for  Rx  parallel  to  X  at  the  0.5  km  separation,  the  processed 
difference  level  is  close  to  the  instrument  noise  level.  The  simple  difference  Ry  -  Y  is  also 
close  to  the  instrument  noise  level,  so  that  the  improvement  due  to  the  processing  is,  on  the 
average,  about  0  dB.  Note  that  Figs.  17  and  18  are  self-consistent;  the  gain  error  in  Rx  ver¬ 
sus  X  produced  a  rather  large  difference  signal  that  was  corrected  by  a  non-unity  value  of 
the  coefficient  A,  so  that  a  clear  processing  improvement  wsb  obtained.  However,  very  little 
gain  error  was  present  between  Ry  and  Y\  hence,  the  coefficient  B  in  Fig.  18  is  near  unity, 
the  simple  first  difference  is  very  close  to  instrument  noise  level,  and  the  processing  yields 
no  improvement. 

Figure  19  gives  results  for  the  14.6  km  separation  in  California  for  data  taken  during 
light  winds.  The  coefficient  A  is  relatively  flat  and  has  a  mean  value  of  0.86;  again  the  evi¬ 
dence  indicates  a  difference  between  the  gains  of  the  X  and  Rx  sensors,  rather  than  an 
alignment  error.  The  B  and  C  coefficients  show  only  a  small  amount  of  structure.  The  large 
variations  in  coefficient  C  are  due  in  part  to  the  fact  that  the  spectral  amplitudes  in  the  Z 
sensor  are  intrinsically  lower  (see  Fig.  9),  so  that  the  "signal”  power  available  to  improve  the 
site-to-site  coherence  is  closer  to  the  sensor  noise  level. 

T  processed  difference  in  Fig.  19  is  approximately  equal  to  the  instrument  noise 
level  be  v-een  2  and  30  Hz  but  exceeds  the  instrument  noise  level  above  and  below  this 
frequency  interval.  As  in  Fig.  17,  the  fact  that  the  processed  difference  cannot  be  sup¬ 
pressed  to  near  the  instrument  noise  level  (as  predicted  by  Eq.  (21))  at  all  frequencies  may 
indicate  the  existence  of  a  factor  that  cannot  be  corrected  by  the  construct  of  Eq.  (1).  The 
improvement  due  to  processing,  also  plotted  in  Fig.  19,  is  on  the  average  about  22  dB  in  the 
2  to  20  Hz  interval.  Most  of  this  improvement  is  probably  due  to  gain  correction  and  correc¬ 
tion  for  any  small  sensor  misalignments,  rather  than  for  any  geology-related  corrections. 

This  plotted  improvement  curve,  we  again  emphasize,  represents  a  comparison  with  the  first 
difference.  A  comparison  between  the  processed  difference  and  the  raw  spectra  (Fig.  9)  indi¬ 
cates  that  the  processing  suppresses  ambient  noise  in  this  frequency  band  by  up  to  40  dB. 

We  believe  the  lack  of  any  structure  or  trend  in  the  coefficients  at  14.6  km  indicates 
that  any  decorrelating  effects  due  to  an  inhomogeneous  geology  are  minimal  at  this  separa- 
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tion.  The  geology  at  China  Lake,  California,  is  relatively  homogeneous— all  Bites  are  within 
the  large  (~650  km2)  Indian  Wells  Valley,  which  has  a  uniformally  layered  level  topology— 
and  the  sites  were  located  at  least  8.0  km  from  the  perturbing  influence  of  the  surrounding 
mountains.  It  appears  then  that  separations  larger  than  14.5  km  are  necessary  to  demon¬ 
strate  the  effects  of  a  nonuniform  geology  at  the  China  Lake,  California  site.  (Data  were 
not  taken  at  larger  separations  because  of  the  limited  transmission  range  of  the  UHF  telem¬ 
etry  link.) 


Figure  20  shows  the  results  of  processing  a  set  of  data  that  produced  definite  structure 
in  the  coefficients,  These  data  were  taken  at  a  separation  of  14.5  km  during  a  period  of 
fairly  high  winds.  The  spectra  and  coherencies  for  this  period  of  time  were  given  in  Fig.  8 
and  14,  respectively.  The  spectra  show  prominent  peaks  between  1  and  5  Hz  that  are  caused 
by  movements  of  the  SQUID  sensors  in  the  earth’s  magnetic  field,  in  spite  of  the  use  of  a 
wind  screen  around  the  instruments,  The  wind  motions  are  uncorrelated,  a  fact  reflected  in 
the  decreased  coherencies  between  1  and  5  Hz  in  Fig.  14.  The  wind  was  gusty,  so  that  decor¬ 
relation  of  the  parallel  sensor  outputs  occurs  only  during  the  gusts;  hence,  the  plotted 
coherence  doeB  not  decrease  entirely  to  zero,  but  is  depressed  by  an  amount  depending  on 
the  relative  duration  of  the  gust  and  quiet  intervals.  The  X  and  Rx  coherence  is  lowered 
more  than  the  Y  and  Ry  coherence.  This  probably  occurred  because  the  wind  was  blowing 
from  the  north,  a  circumstance  that  causes  the  magnetic  north-oriented  X  and  Rx  sensors  to 
produce  a  larger  motion-induced  signal. 

The  motion  noise  in  the  1  to  5  Hz  band  produces  the  structure  in  the  coefficients  in 
Fig.  20.  The  comparison  o t  the  first  difference  and  the  processed  difference  in  Fig.  20  shows 
that  the  processing  improves  the  difference  somewhat  in  the  band  of  motion  noise  inter¬ 
ference,  but  only  to  a  small  degree.  Comparing  Figs.  10  and  20,  the  motion  noise  from  the 
wind  clearly  decreases  the  amount  of  improvement  that  the  processing  can  achieve.  To 
emphasize  this  point,  in  Fig.  21  smoothed  processing  improvement  curves  taken  from 
Figs.  10  and  20  have  been  redrawn  to  show  the  differences.  The  plotted  curves  are  the  ap¬ 
proximate  average  value  of  the  Improvement  with  variability  removed.  Below  about  8  Hz 
the  sample  of  data  contaminated  by  wind  noise  shows  a  marked  decrease  in  processing 
improvement. 


ThlB  observation  of  the  effect  of  motion-induced  noise  suggests  a  source  of  the  uncor- 
rectable  decorrelating  factor  that  causes  the  upturn  in  the  processed  difference  curves  in 
Figs.  17, 18  and  19  at  frequencies  below  2  Hz  (and  that  causes  a  concomitant  drop  in  the 
Improvement-due-to-processing  curves).  A  continuous  background  of  microseismlc  activity 
is  always  present  and  ;•  caused  both  by  natural  worldwide  earthquake,  tidal,  and  ocean  wave 
activity  and  by  local  man-made  activities  (such  as  automobile  and  truck  traffic)  [15] .  This 
microseismic  activity  has  a  power  spectrum  that  increases  in  amplitude  with  decreasing  fre¬ 
quency.  It  is  possible  that  below  a  frequency  of  about  2  Hz  continuous  motion-induced 
noise  caused  by  microseismic  activity  prevents  the  complete  cancellation  of  noise  down  to 
the  SQUID  sensor  noise  level.  We  can  derive  an  approximate  value  for  seismic-induced  mag¬ 
netometer  noise  that  serves  to  demonstrate  at  least  the  possibility  of  this  noise  source  being 
observable  in  the  processing  results.  Iyer  and  Hitchcock  [16]  present  the  results  of  a  seismic 
noise  survey  taken  in  Long  Valley,  California,  which  is  200  km  north  of  China  Lake’s  Indiaik 
Wells  Valley.  The  geological  structure  of  the  two  valleys  is  similar,  consisting  of  a  soft  sedi¬ 
mentary  basin  overlying  a  solid  rock  basement,  so  that  the  seismic  noise  values  in  Long 
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Valley  should  be  Indicative  of  the  seismic  noise  background  at  China  Lake.  In  the  1-  to  2-Hz 
band  the  seismic  background  was  measured  [16]  by  a  vertical  velocity-sensitive  seismometer 
to  vary  from  6  dB  to  24  dB  above  a  reference  level  of  4  X  10"8  m/sA/Hz.  The  range  of 
values  reflects  both  a  spatial  and  temporal  variation  in  noise  level.  We  shall  take  18  dB  (a 
factor  of  63  in  amplitude)  as  a  reasonable  median  value.  Assuming  vertically  polarized 
shear  waves  (vertical  transverse  motion),  or  SV  waves,  the  ground  motion  amplitude 
a  is  related  to  the  rate  A  by  a  ■  £2/2 irf,  where  f  is  the  frequency.  Thus,  at  f  »  1  Hz,  there 
results  a  •  (4  X  10'  8  )(63)/2tr(l)  ■  4  X  10“  7mA/Hi,  During  the  passage  of  an  SV  wave  of 
wavelength  X,  the  peak  tilt  excursion  80  is  related  to  the  wave  amplitude  A  by 


9o 


2ira 
X  ‘ 


(32) 


The  signal  A H  induced  in  the  X  SQUID  sensor  (oriented  magnetic  north-south)  when  ro¬ 
tated  through  a  small  angle  6 0  in  an  earth’s  magnetic  field  H,  that  has  a  dip  angle  </>  is  given 
by  [10] 


A H  -  Htd0  sin  <t> 


2m 


H,  sin  <t> . 


(33) 


At  China  Lake,  <)>  -  60°  and  H.  -  0.B  G.  Values  of  X  were  not  measured  in  the  work  re¬ 
ported  in  Ref.  16,  but  Capon  [17]  quotes  values  of  X  *  4  km  as  typical  for  seismic  SV 
noise  at  1  to  2  Hz.  Substituting  the  values  given  above  into  Eq.  (33)  yields  an  estimate  for 
seismic-induced  noise  of  A H  ■  3  X  10"  10GA/Hz.  This  value  is  of  the  order  of  the  level  of 
residual  noise  at  1  to  2  Hz  in  the  processed  differences  shown  in  Figs.  17, 18,  and  19.  The 
calculation  indicates  that  seismic  noise  may  determine  the  processed  difference  level  in  the 
lower  portion  of  the  0.3  to  40-Hz  band. 

Norway  Data 

Figure  22  displays  the  results  for  the  0,5  km  separation  for  the  Tromso,  Norway  data. 
These  results  were  computed  using  Rx  for  the  remote  site  sensor.  When  the  SQUID  instru¬ 
ment  was  installed  at  the  0.5-km  remote  site,  the  intention  of  course  was  to  place  Rx  paral¬ 
lel  to  X,  Ry  parallel  to  Y,  and  Rt  parallel  to  Z,  However,  we  believe  on  the  basis  of  the  co¬ 
herence  plots  in  Fig,  15  and  the  appearance  of  the  coefficient  plots  in  Fig,  22  that  the  re¬ 
mote  site  was  installed  with  the  horizontal  axes  inadvertently  rotated  by  an  angle  6  from  the 
parallel  alignment.  In  Fig.  22  a  dotted  line  has  been  drawn  at  the  average  value  of  A  and  B. 
For  axes  at  the  remote  site  rotated  in  the  horizontal  plane  by  an  angle  8,  Appendix  B  shows 
that  the  coefficients  A  and  B  are  given  by 
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IAI 


cos  8  -  l7xyl  sin  8 
1-7  2XY 


(34) 


IBI 


sin  8  -  l7Xy-l  cos  6 
1-  7 XY 


(36) 


Appendix  B  emphasizes  the  assumptions  that  enter  into  the  derivation  of  these  equations. 
Given  that  we  are  interested  only  in  an  approximate  value  for  a  rotation  angle,  we  believe 
that  the  assumptions  are  reasonably  well  satisfied.  Figure  15  indicates  that  y\  Y  at  0.08 
(above  5  Hz).  Substituting  A  ■  0.75  and  y xy  m  0.28  into  Eq.  (B-14),  which  is  Eq.  (34) 
solved  for  sin  8,  gives  0  ■  32°.  Similarly,  substituting  B  ■  0.37  and  yxY  m  °-28  in  Eq. 
(B-15),  which  is  Eq.  (35)  solved  for  cos  Q,  gives  $  -  33°.  These  two  values  for  the  angle  of 
rotation  are  in  reasonable  agreement,  in  view  of  the  approximations  used  and  the  uncer¬ 
tainties  in  the  variables.  An  estimate  of  the  angle  8  can  also  be  obtained  directly  from  the 
coherence  plot  in  Fig.  15  by  using  Eq.  (B-l): 


0.60, 


(36) 


where  0.60  is  the  approximate  average  value  of  the  coherence  between  3  and  10  Hz  on 
Fig.  15,  Equation  (36)  gives  a  value  of  8  “  39°.  This  value  is  in  reasonable  agreement 
with  the  values  for  the  angle  derived  from  the  coefficient  values. 

Partially  because  the  average  magnitudes  of  the  coefficients  are  lower,  the  variability 
across  the  frequency  band  is  relatively  large.  Large  excursions  in  the  coefficients  are  also 
noticeable  at  the  frequencies  of  the  tonal  interference  present  in  the  spectra  (Fig.  10). 

The  processed  difference  in  Fig.  22  shows  a  form  similar  to  the  processed  differences 
evaluated  for  the  California  data.  From  2  to  30  Hz  the  difference  is  near  the  SQUID  sensor 
noise  level  but  increases  below  this  interval,  possibly  because  of  microseismic  motion  of  the 
sensor  that  cannot  be  compensated.  The  simple  first  difference  shows  an  unusual  sharp  in¬ 
crease  below  2  Hz,  presumably  related  to  the  30s  misalignment  of  the  sensor  axes.  Because 
of  this  sharp  increase,  the  Improvement  due  to  processing  (relative  to  the  first  difference)  in 
Fig.  22  is  nearly  35  dB  at  0.3  Hz. 

Figure  23  displays  the  processing  results  for  the  9.3  km  separation  in  Norway.  The  co¬ 
efficients  contain  a  large  amount  of  structure  that  was  not  observed  in  the  other  data.  A 
large  peak  in  all  of  the  coefficients  is  noticeable  at  the  frequency  of  the  Schumann  reso¬ 
nance  fundamental  at  about  8  Hz.  We  have  no  explanation  of  this  large  peak  in  the  coeffi¬ 
cients  except  to  note  that  the  spectra  for  this  data  (Fig.  11)  show  an  abnormally  large 
Schumann  resonance  fundamental  peak  (relative  to  the  resonance  harmonics  and  the  over¬ 
all  background  level). 
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Also  of  interest  in  Fig.  23  is  the  large  increase  in  B  at  frequencies  below  6  Hz.  This 
structure  in  the  coefficient  is  probably  related  to  inhomogeneities  in  the  local  geology.  The 
particular  feature  of  the  geology  that  may  cause  this  structure  is  the  conductivity  anomaly 
presented  by  Balsfjord,  shown  in  Fig.  3.  The  form  of  the  coefficient  B  indicates  that  as  the 
frequency  of  interest  is  progressively  decreased  below  7  Hz,  a  larger  portion  of  the  base  site 
Y  axis  (which  is  orthogonal  to  the  Rx  axis  used  in  this  processing)  must  be  admixed  to  the 
difference  signal  formed  according  to  Eq.  (1).  The  conductivity  anomaly  presented  by  the 
fjord  thus  appears  to  rotate  the  incident  noise  fields  at  the  remote  site  relative  to  the  base 
site. 


The  processed  difference  in  Fig.  23  falls  at  approximately  the  instrument  noise  level 
across  the  entire  frequency  Interval.  The  increase  in  the  difference  level  at  frequencies  below 
2  Hz  is  not  as  pronounced  as  the  increase  in  the  processed  differences  given  in  Figs.  17 
through  20,  and  in  Fig.  22.  The  improvement  in  the  difference  due  to  the  processing  has  a 
maximum  value  of  about  20  dB  at  10  Hz.  Comparing  the  raw  spectrum  in  Fig.  11  for  the 
X  base  sensor  with  the  processed  difference  indicates  that  ambient  noise  is  suppressed  by 
35  dB  by  the  processing. 


CONCLUSIONS  AND  RECOMMENDATIONS 

No  substantial  differences  were  observed  between  the  California  data  and  the  Norway 
data  in  either  the  characteristics  of  the  noise  or  the  ability  of  the  FSB  algorithm  to  suppress 
the  noise.  In  both  locations  the  FSB  algorithm  was  able  to  reduce  the  ambieut  noise  essen¬ 
tially  to  the  SQUID  instrumental  noise  level.  A  small  increase  in  the  processed  difference 
noise  level  below  5  Hz  is  probably  a  result  of  seismic-induced  motion  noise  for  which  the 
FSB  processing  cannot  compensate. 

For  a  detection  application  of  the  FSB  processing,  the  reduction  in  ambient  noise 
achievable  when  the  processed  difference  is  evaluated  using  fresh  time  series  data  with  coef¬ 
ficients  determined  many  hours  earlier  must  be  investigated.  We  have  not  yet  processed  the 
data  in  this  manner  but  hope  to  address  this  point  in  a  future  report. 

The  usefulness  of  component  differencing  as  a  geophysical  exploration  technique 
depends  on  the  ease  with  which  the  specific  form  of  the  coefficients  can  be  interpreted 
(with  the  aid  of  a  mathematical  model)  in  terms  of  subsurface  geological  features.  The 
development  of  the  theory  to  explain  the  specific  form  of  the  coefficients  (and  the  coher¬ 
ence)  measured  at  the  9.3-km  separation  in  Norway,  where  geological  effects  due  to  the 
fjord  clearly  were  observed,  is  a  worthwhile  area  for  further  work. 

Finally,  the  numerous  adaptive  noise  suppression  techniques  described  in  the  literature 
by  Wldrow  [IS] ,  Griffiths  [19] ,  and  others  should  also  be  applied  to  the  data  taken  in  this 
study  for  comparison  with  the  FSB  algorithm.  We  hope  to  make  this  comparinon  ns  a  con¬ 
tinuation  of  the  present  study. 
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Appendix  A 

PROCESSED  DIFFERENCE  NOISE  LEVEL 


This  appendix  derives  an  expression  for  the  noise  level  of  the  processed  difference 
signal  A.  Explicit  computation  of  the  noise  level  of  A  is  very  complicated  for  the  general 
3-dimensional  case,  and  hence  we  calculate  it  for  a  2-dimensional  case  ( Z  ■  0).  The 
2-dimensional  calculation  is  sufficient  for  deriving  an  expression  to  compare  with 
experiment. 

The  equation  that  A  and  B  should  satisfy  is  of  course  given  by 

A^R-  AX- BY,  (A-l) 

We  now  seek  the  average  noise  power  in  A,  denoted  o?,  when  A  and  B  (as  evaluated  in 
Eqs.  (8)  and  (9))  are  used  with  a  series  of  measurements  of  R,  X,  and  F: 


N 


N 


«I  ■  -  AXI  -  *  W*  -  A‘x;  -  B-  y* ). 

I  i 

Substituting  Eqs.  (8)  and  (9)  into  Eq.  ( A-2)  and  making  use  of  the  identity 

V  +  t)  if  +  8*  )m  l/l 2  +  Itfl 2  +  2  Re  { fg }  , 
the  following  expression  results  after  tedious  algebra: 

1-  Irx"  Tly"7xy  +  2Re{TflxTfly7Xy}' 


a\  -  o%  I  - 


(A-2) 


(A-3) 


(A-4) 


where  the  notation  yfJ  ■  coh  (IJ)  (/  -  R,  X\  Jm  X,Y)  has  been  used.  Several  '.nlting  forms 
of  Eq.  (A-4),  and  their  implications,  are  examined  in  the  text  under  Expected  Noite  Level  of 
the  Difference  Signal. 

We  now  develop  an  expression  for  the  coherence  that  separates  the  effect  of  the 
uncorrelated  intrinsic  sensor  noise  from  the  portion  of  the  signals  measured  at  each  site 
that  are  possibly  coherent.  This  expression  is  to  be  substituted  into  Eq.  (A-4)  to  obtain 
an  estimate  of  c^.  We  assume  that  the  measured  FT’s  can  be  written  as: 


X{f)-Xc(f)+Xn(f), 

Y{f)  -  Yc(f)  +  Yn(f),  and 

m-Rc(f)  +  Rn(f), 
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(A-5) 


t., 
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where  the  subscript  c  refers  to  a  “true”  signal  that  may  be  correlated  at  the  base  and 
reference  sites,  and  the  subscript  n  refers  to  uncorrelated  noise  components  arising,  for 
example,  from  the  intrinsic  sensor  noise.  If  we  use  as  an  example  X  and  R  from  Bq.  ( A-5), 
the  coherence  function  7^  of  Eq.  (26)  can  be  written  (dropping  the  explicit  frequency 
dependence) as 


where 


~2 


XCRC 


(*c*c>  <Mc> 


(A-6) 


(A-7) 


is  the  true  coherence  between  the  signals.  As  one  would  expect,  the  additive  sensor  noise 
causes  the  measured  coherence  to  be  less  than  the  true  coherence.  Similar  equations  can  be 
written  for  7y^  and  7^y.  The  ratios  and  (o%n!OftC)iri  the  denominator  of  Bq. 

(A-6)  are  the  noise  power  to  signal  power  ratios  of  the  sensors  (signal  here  refers  to  the  true 
external  signals  Incident  on  the  SQUID  magnetometers);  as  will  be  seen  in  the  data  discussed 
below,  these  ratios  are  approximately  equal  for  all  sensors  and  are  no  larger  than  about 
0.01.  Therefore,  let  these  ratios  be  denoted  6.  Then  Eq.  (A-6)  can  be  written: 


1  +  25+5® 


‘^(1-26), 


(A-8) 


where  the  fact  that  6  1  has  been  used  to  obtain  the  approximation.  If  Eq.  (A-8)  and  a 

similar  equation  for  7 are  substituted  into  Eq.  (A-4),  the  following  expression  results, 
where,  for  simplicity,  72  y.  -  0  has  been  assumed: 

°A  "  °fi  ~  7r  Y  +  26  (  Tj?x  +  y  ) ]  •  A'9) 


Again  viewing  the  noise  cancellation  process  as  the  measurement  of  an  input  (X  and  Y )  - 
output  (R)  linear  system  transfer  function,  [8]  shows  that  if  the  system  is  strictly  linear, 
then: 


yRx  +  yRY 


(A-10) 
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always  This  suggests  that  the  existence  of  any  nonlinear  effects  und  unconeiated  noise 
sources  other  than  the  intrinsic  sensor  noise  can  quantified  by  a  parameter  e  by  letting: 


e  m  1 "  x  ~1ry  •  (A'H) 

Substituting  Eq.  (A-ll)  into  Eq.  (A-9)  gives 

o£-  c/2  (e  +  26),  (A-12) 


where  the  term  2Se  has  been  dropped  on  the  assumption  that  it  is  a  second  order  term. 
Noting  that 


JRn 


JRn 


aRC  +  °Rn 


(A-13) 


Eq.  (A-12)  can  be  written 


1  2  a 


Rn 


+  e  oi 


(A-14) 


Equation  (A-14)  is  the  simplest  possible  equation  for  the  noise  power  in  the  difference 
signal  and  demonstrates  that ’the  FSB  differencing  technique  decreases  the  noise  power  to 
approximately  the  noise  level  of  the  sensors  (the  factor  of  2  arises  from  the  fact  that  the 
sensors  are  used  pair  wise)  plus  a  residual  amount  arising  from  nonlinearities  or  uncorrelated 
noise  sources.  Although  y%y  was  set  equal  to  zero  to  go  from  Eqs.  (A-4)  to  Eq.  (A-14),  a 
finite  value  of  y%  y  makes  Eq.  (A-14)  somewhat  more  complicated  but  doeB  not  change  the 
conclusions. 


Appendix  B 

THE  EFFECT  OF  AN  ANOULAR  MISALIGNMENT  ERROR 

This  appendix  derives  an  expression  for  the  coefficients  A  and  B  when  the  base  site 
ho  rizontal  axes  and  the  remote  site  horizontal  axes  are  rotated  with  respect  to  each  other  by 
an  angle  $.  The  expressions  are  derived  by  making  :  he  following  simplifying  assumptions: 

•  Only  the  horizontal  axes  are  considered,  i.e.,  a  two-dimensional  case  with  Z*  0 
is  used. 

*  Only  decorrelating  effects  from  misalignment  «re  considered,  l.e,,  geological  inhomo¬ 
geneities  are  nut  explicitly  included  in  any  manner. 

Equations  (8)  and  (9)  given  an  explicit  expression  for  A  and  B  in  the  two-dimen¬ 
sional  cast.  The  misalignment  of  the  axes  at  the  base  and  remote  sites  causes  the  coherence 
functions  that  appear  in  these  equations  to  be  reduced.  We  assert  that  the  dependence  of  the 
coherence  on  a  presumed  misalignment  0  between  the  X  and  R  axes  is  given  approximately 
by 


Ieoh  (XR)1 »  cos  6  . 


(B-l) 


We  can  demonstrate  this  dependence  heuristlcally  in  the  following  manner.  The  cross¬ 
spectra  that  enter  into  the  numerator  of  the  coherence  function  defined  in  Eq.  (10)  or  Eq. 
(26)  can  be  represented  by  [8] 


xm<f) 


lim  lim  ^ 
Af-*o  T-*<*>(Af)T 


T 

f  X{t,fAf)R*(t,fAf)dt, 

C 


(B-2) 


where  X{t,f,i\f)  and  R(t,f,Af)  are  X  and  R  time  series  that  have  been  filtered  through  a 
(narrow)  band  pass  filter  A/,  The  time  series  X(t)  and  R(t)  for  ULF  and  ELF  noise  can  be 
represented  as  the  superposition  of  many  overlapping  random  wave  packets,  as  shown  in 
Fig.  B-l;  these  packets  arise  from  the  filtering  at  a  relatively  low  frequency  f  of  atmospheric 
impulses  whose  inverse  pulse  width  is  much  larger  than  f.  R  eferring  to  Fig.  B-2,  we  can  then 
write  the  signals  R(t)  and  X{t)  as  follows: 


X(t)  -  J]  F( T,,Z>1.,tf|.)e'2ff't'  cos  ^  , 

i 


(B-3) 
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Pig.  B-l  —  Generalized  representation  of  ELF  wavepacketi  that  aum  togathar  to 
form  tha  recorded  notae  voltage 
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Fig,  B-2  —  Geometry  of  X  and  R  aenaora  with  a 
mlaallgnment  of  6,  and  plane-polarized  Incident 
wave  at  angle  \f 


m  -  £  F(TitDitHt)J**f* cos  (*,  -  6) , 


where  F(ri,Dj,Hj)  is  a  generalized  representation  of  the  individual  wavepacket  envelopes 
with  start  time  rh  duration  Dh  and  peak  amplitude  Ht.  As  Fig,  B-2  indicates,  ^  is  the  azi- 
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muthal  arrival  angle  of  the  packet.  The  quantities  r,,  Dit  Hjt  and  ^  are  assumed  to  be  ran¬ 
dom  variables 

Forming  the  product  needed  in  Eq.  (B-2)  gives 


XR  “  ]££ F(ritDitHi)F\rk,DhtHk)  cos  p,  cou(*k  -  6) .  (B-B) 

I  k 


The  cross  terms  in  the  summation  will,  on  the  average,  sum  to  zero  because  of  the  ran¬ 
dom  nature  of  the  wave  packets.  There  remains: 


XR  -  cos  ft  cos  (ty  -  0) . 

i 


(B-6) 


The  denominator  in  Eq,  (26)  for  the  coherence  is  found  in  a  similar  manner.  The  two 
terms  are: 


XX  -  COS2  ifi(  , 

t 


(B-7) 


and 


RR  -  cos2  ty  -  0) . 

i 


(B-8) 


The  computation  of  the  coherence  requires  the  averaging  of  sequential  terms  of  the  form  of 
Eqs.  (B-6),  (B-7),  and  (B-8).  The  average  is  taken  over  all  of  the  random  variables  in  the 
terms.  Thus,  the  averaging  amounts  to  converting  the  summation  to  an  integral  over  rr  D., 
Hh  and  The  equation  for  the  coherence  can  be  written 


Ieoh  (Xft)|  - 


2n 


cos  i p  cos(^>  -  0  )d>fi 


2rr 

1/2 

"  2tt 

fF2drdDdH  f  cos 2<fid<p 
J  •/ 

jF2drdDdH  f  cos 2(v?-0)dy> 

O 

O 

(B-9) 
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The  integrations  over  t,  D,  and  H  produce  terms  that  cancel  in  the  numerator  and  denom¬ 
inator;  there  remains 


Ieoh  (ATJ?)| 


2rr 

/« 


cos  cob  6)difi 


2rr 

/» 

1/2 

2v 

/  cos2 

1  cos2  (t fi-9)d\fi 

J 

L  0 

(B-10) 


Performing  the  integrations  yields 


|coh  (XR) |  “  cos  6  , 


which  is  Eq.  (B-l). 

Since  we  are  assuming  that  there  are  no  decorrelating  effects  due  to  geological  inhomo¬ 
geneities  (and,  by  implication,  that  sensor  noise  is  negligible),  Eq.  (A-10)  implies  that 


Ieoh  (YB)|  -  sin  9  , 


(B-H) 


Then,  substituting  Eqs.  (B-l),  (B-ll),  and  the  definition  I7*  y|  -  |coh(Xy)|  into  Eqs.  (8) 
and  (9)  yields 


A  » 


cos  9  -  i7xyl  sin  9 

1  -yxv 


and 


Ifll 


sin  9  -  \yXy\  cos  9 


(B-12) 


(B-l  3) 
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These  are  the  desired  equations.  In  this  report  values  of  A  and  B  are  determined  from 
the  recorded  data  and  then  used  to  calculate  6.  For  this  purpose,  fiqs.  (B-12)  and  (B-13)  are 
solved  for  sin  8  and  cos  8  respectively: 

„  ^xy ~  ^xya  ±  jjMryd “  7*y) +  ^xy^xy  +^2)"  A2  + 1]  12 
•In  8  - - - - - - - —  ,  (B-14) 

1  *^XY 


and 


.  <2 XY  ~  1)7 xyB  ±  [27*yd  -  7&y)  +  7 %y(V$y  *  Ba)  -  £a  +  lj12 
cos  8  m  ————— — — — — . - — 1  -  —  „  —  —  

1+7  $y 


(B-16) 
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